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We construct the operator that projects on the physical states in loop quantum gravity. To 
this aim, we consider a difTeomorphism invariant functional integral over scalar functions. The 
construction defines a covariant, Feynman-like, spacetime formalism for quantum gravity and relates 

OO ' this theory to the spin foam models. We also discuss how expectation values of physical quantity 

0^ , can be computed. 
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■ The loop approach to quantum gravity ^ based on the Ashtekar variables |^ has been successful in establishing 
\ a consistent and physically reasonable framework for the mathematical description of quantum spacetime ||^. This 
. framework has provided intriguing results on the quantum properties of space, most notably detailed quantitative 

■ results on the discrete quanta of the geometry jj]. The nonperturbative dynamics of the quantum gravitational 
^ field, however, is not yet well understood. Two major questions are open. First, several versions of the hamiltonian 

I constraint have been proposed but the physical correctness of these versions has been questioned Sec- 

Cn| . ond, a general scheme for extracting physical consequences from a given hamiltonian constraint, and for computing 
expectations values of physical observables is not available. 

In this paper, we address the second of these problems - a solution of which, we think, is likely to be a prerequisite 
for addressing the first problem (the choice of the correct hamiltonian constraint). The problem we address is how 
expectation values of physical observables can be computed, given a hamiltonian constraint operator H{x). (For an 
earlier attempt in this direction, see |^.) We address this problem by constructing the "projector" P on the physical 
Hilbert space of the theory, namely on the space of the solutions of the hamiltonian constraint equation. Formally, 
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I this projector can be written as 

bj:). 



P ^l[SiHix)) ^ J [DN] e"^/'*'" m^Wi^) (i) 



^ . in analogy with the representation of the delta function as the integral of an exponential. The idea treating first class 
■ constraints in the quantum theory by using a projector operator defined by a functional integration has been studied 
by Klauder , and consider also by Govaerts , Shabanov and Prokhorov , Henneaux and Teitelboim jl^] and 
others. In the case of gravity, the matrix elements of P between two states concentrated on two 3-geometries g and 
g' can be loosely identified with Hawking's propagator P{g,g'), which is formally written in terms of a functional 
integral over 4-geometries ||TJ] . 

A step towards the definition of the projector P was taken in where a perturbative expression for the exponential 
of the hamiltonian smeared with a constant function N(x) = N was constructed. What was still missing was a suitable 
diffeomorphism invariant notion of functional integration over N{x). Here, we consider an integration on the space 
of the scalar functions N{x). This integral, modeled on the Ashtckar-Lewandowski construction ||l5| and considered 
by Thiemann in the context of the general covariant quantization of Higgs fields pGf , allows us to give a meaning to 
the r.h.s. of (|^). Using it, we succeed in expressing the (regularized) matrix elements of the projector P in a well 
defined power expansion. We then give a preliminary discussion of the expectation values of physical observables. 
The construction works for a rather generic form of the hamiltonian constraint, which includes, as far a we know, the 
various hamiltonians proposed so far. 

As realized in , the terms of the expansion of P are naturally organized in terms of a four-dimensional Feynman- 
graph-like graphic representation. The expression (^ can thus be seen also as the starting point for a spacetime 
representation of quantum gravity. Here, we complete the translation of canonical loop quantum gravity into covariant 
spacetime form initiated in The "quantum gravity Feynman graphs" are two-dimensional colored branched 
surfaces, and the theory takes the form of a "spin foam model" in the sense of Baez or a "worldsheet theory" in 
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the sense of Reisenberger , or a "theory of surfaces" in the sense of Iwasaki , and turns out to be remarkably 
similar to the Barret-Crane model and to the Reisenberger model (see also [^). On the one hand, the 
construction presented here provides a more solid physical grounding for these models; on the other hand, it allows 
us to reinterpret these models as proposals for the hamiltonian constraint in quantum gravity, thus connecting two of 
the most promising directions of investigations of quantum spacetime pst - 

The paper is organized as follows. In Section ||, the basics of loop quantum gravity are reviewed, organized from 
a novel and simpler perspective, which docs not require the cumbersome introduction of generalized connections, or 
projective limits (see also p^ ). Section [II presents the definition of the diffeomorphism invariant functional integral. 
In section we construct the projector P and discuss the construction of the expectation values. 



II. LOOP QUANTUM GRAVITY 

General relativity can be expressed in canonical form in term of a (real) SU (2) connection A defined over a 3d 
manifold E H,^- We take S to be topologically 5*3. The dynamics is specified by the usual Yang-Mills constraint, 
which generates local SU{2) transformation, the diffeomorphism constraint D[f], which generates diffeomorphisms 
of S, and the Hamiltonian constraint H[N], which generates the evolution of the initial data in the (physically 
unobservable) coordinate time. Here / is in the algebra of the group Diff of the diffeomorphisms of E, namely it is a 
smooth vector field on E, and iV is a smooth scalar function on E. The theory admits a nonperturbative quantization 
as follows. (For a simple introduction, see for details see ^ and references therein.) 



A. Hilbert space and spin networks basis 

We start from the linear space L of quantum states ^'(A) which are continuous (in the sup-topology) functions of 
(smooth) connections A. A dense (in L's pointwise topology) subset of states in L is formed by the graph-cylindrical 
states A graph-cylindrical state ^r,/(^) is a function of the connection of the form 

^rjiA)^f{U{ei,A),...,UienA)), (2) 

where F is a graph embedded in E, ei . . . e„ are the links of F, U (e, A) ~ Vexp A is the parallel propagator SU (2) 
matrix of A along the path e, and / is a complex valued (Haar-integrable) function on [SU{2)]"'. The function '^rj{A) 
has domain of dependence on the graph F; one can always replace F with a larger graph F' such that F is a subgraph of 
F', by simply taking / independent from the group elements corresponding to the links in F' but not in F. Therefore 
any two given graph-cylindrical functions can always be viewed as defined on the same graph F. Using this, a scalar 
product is defined on any two cylindrical functions by 

(^'r,/,*r.3) = / dUi... dUn f{Ui . . . U„) g{Ui . . . Un) (3) 

J[S(7(2)]" 

and extends by linearity and continuity to a well defined scalar product on L. The Hilbert completion of L 

in this scalar product is the Hilbert space Haux- the quantum state space on which quantum gravity is defined]^ We 
refer to p4| for the construction of the elementary quantum field operators on this space. 

The SUl2) gauge invariant states form a liner subspace Lq in L. A convenient orthonormal basis in Lq is the spin 
network basis [Q, constructed as follow. Consider a graph F embedded in E. To each link e of F, assign a nontrivial 
SU{2) irreducible representation je, which we denote the color of the link. Consider a node n of F, where the links 
ei . . . ejv meet; consider the invariant tensors v on the tensor product of the representation je^ . . . je„ of the links that 
meet at the node; the space of these tensors is finite dimensional (or zero dimensional) and carries an invariant inner 
product. Choose an orthogonal basis in this spac^, and assign to each node n of F one element w„ of this basis. 



^Haux is the state space of the old loop representation equipped with a scalar product which was first obtained through a 
path involving C*-algebraic techniques, generalized connections and functional measures 10261 . Later, the same scalar product 
was defined algebraically in directly from the old loop representation. The construction of Haux given here is related to 
the one in [ p^Jj^ ] but does not require generalized connections, infinite dimensional measures or the other fancy mathematical 
tools that were employed at first. 

■^For later convenience, we choose a basis that diagonalizes the volume operator M,p4|. 
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A spin network S = (F, {j}, {v}) is given by a graph T and an assignment of a color je to each Unk e and a basis 
invariant tensor w„ to each node n. 

The spin network state 4's(A) is defined as 

'^siA)=l[l[vnR'^{Uie,A)) (4) 

e n 

where W {U) is the matrix representing the SU{2) group element U in the spin-j irreducible representation, and 
the two matrix indices of R^^{U{e,A)) are contracted into the two tensors w„ of the two nodes adjacent to e. An 
easy computation shows that (with an appropriate normalization of the basis states w„ p7| ) the states form an 
orthonormal basis in Haux 

= '5r,r'%}{/}'5{t,}{t,'}- (5) 



B. Diffeomorphisms 

The Hilbert space Haux carries a natural unitary representation U{Dif f) of the diffeomorphism group of S. 

cl>eDtff. (6) 

It is precisely the fact that Haux carries this representation which makes it of crucial interest for quantum gravity. 
In other words, Haux and its elementary quantum operators represent a solution of the problem of constructing a 
representation of the semidirect product of a Poisson algebra of observables with the diffeomorphisms. |2g| ] 
Notice that U sends a state of the spin network basis into another basis state 

[Ui(t>)iPs]{A) = iPsi^-'A) = ^^s{A). (7) 

Intuitively, the space H^i // of the solutions of the quantum gravity diffeomorphism constraint is formed by the states 
invariant under U. However, no finite norm state is invariant under U, and generalized-state techniques are needed. 
We sketch here the construction of Hdi / / , because the solution of the hamiltonian constraint will be given 

below along similar lines. Hdiff is defined first as a linear subset of L*, the topological dual of L. It is then promoted 
to a Hilbert space by defining a suitable scalar product over it. Hdiff is the linear subset of L* formed by the linear 
functionals p such that 

pmm = pw (8) 

for any cj) e Diff. From now on wc adopt a bra/ket notation. We write (||) as 

= (pIV) (9) 

and we write the spin network state ^'s as 15). 

Equivalence classes of embedded spin networks under the action of Diff are denoted as s and called s-knots, or 
simply spin networks. We denote as s{S) the equivalence class to which S belong. Every s-knot s defines an element 
(s\ of Hdiff via 

{s\S) =0, if S7^s(S) 

= c^, if s = s(S). (10) 

Here Cs is the integer number of isomorphisms (including the identity) of the (abstract) colored graph of s into itself 
that preserve the coloring and can be obtained from a diffeomorphism of E. A scalar product is then naturally defined 
in Hdiff by 

{s\s') ^ {s\S') (11) 

for an arbitrary S' such that s{S') — s' . One sees immediately that the normalized states -^\s) form an orthonormal 
basis. 

The space Hdiff is not a subspace of Haux (because diff invariant states have "infinite norm"). Nevertheless, an 
important observation is that there is a natural "projector" H from Haux to Hdiff 
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U:\S)^\s{S)) (12) 

which sends the state in Haux associated to an embedded spin network S into the state in Hdiff associated to the 
corresponding abstract spin network state s. Notice that 11 is not reaUy a projector, since Hdiff is not a subspace of 
Haux, but we use the expression "projector" nevertheless, because of its physical transparency. Since Hdiff can be 
seen as a subspace of H*^^,^., the operator 11 defines a (degenerate) quadratic form ( | )diff on Haux 

{S\S')d^ff ^ {S\U\S') = {s{S)\S') = {s{SMS')). (13) 

Hdiff is can be defined also by starting with the pre-Hilbert space Haux equipped with the degenerate a quadratic 
form ( I )diff, and factoring and completing the in the Hilbert norm defined by { \ )diff ■ |@ That is, states in Haux 
are the (limits of sequences of) equivalence classes of states in Haux under { \ )diff- In other words, knowing the 
"matrix elements" 

{S\U\S') (14) 

of the projector 11 is equivalent to having solved the diffeomorphism constraint. 

Furthermore, the above construction can be expressed also in terms of certain formal expressions, which are of 
particular interest because they can guide us in solving the hamiltonian constraint. Define a formal integration over 
the diffeomorphism group Dif f satisfying the two properties 

D^^l, (15) 

Diff 

and 

/ D(j) 5s,ci,s = Cs(s)- (16) 

JDiff 

Then a diff invariant state | s) can be written as a "state in Haux integrated over the diffeomorphism group" . That is 



\s{S)) = / Del) \U{<P)S), (17) 

JDiff 

In fact, the equations ( |lo| ) and (jll|) can be obtained from the equations (|l^), ( [l6| ) and (|l^). Using this, we can write 
the projection operator 11 , defined in (p^), as 

n= /" D(l)U[(j)]. (18) 

Equivalently, we may write the group element as an exponential of an algebra element, and formally integrate over 
the algebra rather than over the group, that is 



U = jDfe-'^^fK (19) 

This equation has a compelling interpretation as the definition of the projector on the kernel of the diffeomorphism 
constraint operator Da{x) via 

U^Yl 5{Da{x)) ^ J Dfe-'I'''^ r(^)Da(x) (20) 



as m 



Six) = ^ 1 dp e-'f^ (21) 
We shall define the projector on the kernel of the hamiltonian constraint in a similar manner. 
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III. A DIFFEOMORPHISM INVARIANT MEASURE 



In this section, we construct a measure on the space of scalar functions [[l6| , which will be needed for defining the 
analog of Eq. ^ for the hamiltonian constraint. Consider smooth functions N : ^ Si on the three- manifold S, 
taking value on the interval / = [0, T[. We keep track of the "length of the interval /" , T, instead of normalizing it 
to one, because this will simplify keeping track of dimensions in the physical application. Let TV be the space of such 
functions, equipped with the sup topology. Let F(N) be a continuous complex function on the infinite dimensional 
topological vector space J\f, and denote the space of these functions as C. Let {xi} = xi, . . . , a;„ be a set of (disjoint) 
points in S, and / : (/)" ^ C a complex integrable function of n real variables. Consider a function F e £ of the 
form 

F^,^}j{N) = f{N{xi),...,N{x.n)), (22) 

namely a function of N having the the set {xi} as its domain of dependence. The set of functions of this form form 
a dense linear subspace of C, in the pointwise topology. 

The simplest nontrivial of such functions is obtained by picking a single point x and choosing f{N) — N . Notice 
that this defines precisely the Gel'fand transform F.j.,{N) = N{x), or in Gel'fand's enchanting notation 

x{N) = Nix). (23) 

Since the functions of the form (|2^) can be seen as a generalization of Gel'fand's Fx{N),we denote them as "generalized 
Gel'fand functions", or simply Gel'fand functions. Gel'fand functions can be seen as the scalar-field analog of the 
Ashtekar-Lewandowski's graph-cylindrical functions (H), which are defined for connection-fields. 
Define the following linear form on the Gel'fand functions 

DN i^{,.},/(iV) = j^^ dm... dN,, f{Ni, Nr.). (24) 

Here is the normalized invariant measure on the interval /. Finally, denote the closure of C in the norm 

\\F\\^ J DN\F{N)\ (25) 

as Li[J\f]; the linear form ( p^ ) extends by continuity (in the topology defined by this norm) to all of Li[J\f]. 
A simple class of integrable functions is given by polynomial Gel'fand functions. We have indeed 

DN l^^ J^dN 1 = 1, 

DN N{x) = ^ J^dN N = ^^T' = (26) 

Notice that for quadratic functionals we must distinguish two cases 

DN N{x) N{y) = ^ JJ^ dN^dN^ iVi iVa = ^ T^, (27) 

DN N{x)N{x) = ^ jdN N^ = ^^T' = ^T\ (28) 

Namely, we must distinguish the case in which the arguments of the two functions N{ ) are distinct or the same. 

ffc) (k) 

The general pattern should be clear. A general polynomial functional Fn-^,,,nK ^il^ have points x\ . . . Xn^ in its 
domain of dependence in which the function N(x) appear with power k. A simple calculation yields then 

^n\...nK — j Fn\...nK 

= JdN [«(4'>)...lV(<))][iV(4^>)...JV(ig)l^..|JV(iS'f'). ..«(<>)]'' 
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The difFeomorphism group Dif f of S acts naturally on A^, via {(j)N){x) = N{(j) ^{x)), where : E — > E is in Diff. 
It is easy to see that the integral is diffeomorphism invariant 

J [DN] F[(I)N] = J [DN] F[N]. (30) 

This follows from the fact that the r.h.s. of ( p^ is clearly insensitive to a diffeomorphism transformation on N. 

IV. DYNAMICS: THE REGULARIZED PROPAGATOR 

We now come to the construction of the physical state space Hphys and the partition function of the theory. We 
have to solve the Dirac's hamiltonian constraint equation 

H[N]^l^ = (31) 

for the quantum Hamiltonian constraint H[N]. 

A. The Hamiltonian constraint: first version 

The operator H[N] that we consider is a small modification of the Riemanian hamiltonian constraint defined in 
[^. However, we make only use of the general structure of this operator, which is common to several of the proposed 
variants. We take a symmetric version of H[N], which "creates" as well as "destroying" links. The matrix elements 
of H\N] are given by 

{MH[Nm = {^\C[N]\<t>) + M^pVpy (32) 

where C[N] is the non-symmetric Thiemann's constraint. 

We recall that the operator C[7V], acting on a spin network state |s), is given by a sum of terms, one per each node 
i of s. Sketchy (a more precise definition will be given below), each such term creates an extra link Cadded joining two 
points, i' and i", on two distinct links adjacent to i, and alters the colors of the links between i and i' and between i 
and i" . The result is multiplied by a coefficient depending only on the colors of s, and by the value of the smearing 
function TV "in the point where the node i is located". This is illustrated in Figure 1. 




FIG. 1. Action of the hamiltonian constraint on a trivalent node. 

It is important to observe that C[N] is defined as a map from Hdiff to H*^^. In this definition there is a subtle 
interplay between diff- invariant and non-diff-invariant aspects of the hamiltonian constraint, which is a key aspect of 
the issue we are considering, and must be dealt with with care. The reason C[N\ is defined on Hdiff, namely on the 
diffeomorphism invariant states, is that it is on these states that the "precise position" of the points i and i" and of 
the link Cadded is irrelevant]^ However, C[N] is not diffeomorphism invariant, and therefore C[iV]|s) is not in Hdiff, 
because a diffeomorphism modifies N . The feature of C[A^]|s) that breaks diffeomorphism invariance is the fact that 
it contains a factor given by the value of Nix) in the point in which the node i is located: this location is not a 
diffeomorphism invariant notion. 

Before presenting the precise definition of C[N], which takes care both of its diff-invariant and its non-diff-invariant 
features, we need to define certain peculiar elements of H*^^, which will appear in the definition. Consider an s-knot 



^More precisely, only on these states can the regulator used in the quantization of the classical quantity be removed. 
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s and let i be one of its nodes. Let iV be a scalar function on E. We define the state |s; i, N) in -ff*„j, by (again, we 
interchange freely bra and ket notation) 

{s;i,N\S) = N{xs,^) {s\S). (33) 

where xsa is the position of the node of S that gets identified with the node i of s in the scalar product. Notice that 
the state \s;i, N) is "almost" a diff-invariant s-knot state: in facts, it is "almost" insensible to the location of S. The 
only aspect of this location to which it is sensible is the location of the node. In fact, in the r.h.s. of (^3|), the diff 
invariant quantity {s\S) is multiplied by the value of N( ) in the point in which the node i is located. The Hamiltonian 



constraint defined by Thiemann acts on diff-invariant states and creates states in of the form (33) 

The precise definition of C[N] is the following 

C[N]\s) = Afis) |sf;^,7V). (34) 

Let us explain our notation. Sum over the repeated indices (3 and i is understood. The index i runs over the nodes 
in s. The index /3 = {I' ,1" ,e' ,e") runs over the couples {I', I") of (distinct) links adjacent to each node i and over e' 
and e", which take the values +1 or —1. The s-knot sf was introduced in |^^. It is defined as the right hand part of 
Figure 1. That is, by adding two new nodes i' and i" on the two links I' and I" (determined by /?) respectively, adding 
a new link Cadded colored 1/2 joining i' and i" , and altering the color of the links joining i' and i (and, respectively, 
i" and i) by +1 or —1 according to the value of e' (respectively e"). A^{s) is a coefficient defined in ||] whose explicit 
form is computed in [^l|] . 
Summarizing, we have 

{S\C[N]\s) = A^{s)Nixs,) {S\sf). (35) 

(where {S\s) = (sis') {s\s{S)}; see and (|l|).) Clearly, C[N] can equivalently be viewed as an operator from 
Haux to -ffdi//, by writing 

{s\C[N]\S)=A^is)Nixs^O (s|s(5)f). (36) 
(Recall that s{S) is the s-knot to which S belongs.) This fact allows us to define the symmetrized operator ( |3^ by 

{S\H[N]\s}^A^is)Nixs,.) {S\s^)+A^{s)Nixs,.) {siS)^\s). (37) 

Only one of the two terms in the r.h.s. of this equation may be non-vanishing: the first, if S has two nodes more than 
s; the second, if S has two nodes less than s. 

We can simplify our notation by introducing an index a — (/?, ±1) where -1-1 indicates that a link is added and —1 
indicates that the link is removed. We obtain 

{S\H[N]\.s) = Alis) N{xs,d {S\s?) (38) 

where 

A^-\s)^A^-^\sf). (39) 

Notice that the precise position of Cadded (and thus i' and i") drops out from the final formula, because of the 
diff invariance of the quantity ('S'|s").Q This is essential, because if a specific position for Cadded had to be chosen, 
diffcomorphismiiiiijw invariance would be badly broken. One can view the coordinate distance e between i and i' (and 



between i and i") as a regulating parameter to be taken to zero after the matrix elements (36) have been evaluated. 
The limit e — > is discontinuous, but {s\Sf) is independent from e for e sufficiently small, and therefore the limit of 
these matrix elements is trivial. Thus, the operator is defined thanks to two key tricks: 

1. the diff invariance of the state \s) acted upon allows us to get rid of the precise position of Cadded] 

2. the lack of diff invariance of (5*1, allows us to give meaning to the point X5 j "where the node is located", and 
therefore allows us to give meaning to the smearing of the operator with a given (non diff invariant!) function 
N{x). 

This is why H[N] is defined as a map from Hdiff to H*^^. The first of these two key facts, which allow the quantum 
hamiltonian constraint operator to exist, was recognized in |p3|, the second in pi. 



*One must only worry about the positioning of Eadded up to isotopy. This is carefully defined in [g], following a construction 
by Lewandowski. 
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B. The Hamiltonian constraint: second version 



The interplay between difF invariant and non difF invariant constructs described above needs to be crafted even 
more finely, in order to be able to exponentiate the Hamiltonian constraint and derive its kernel. In fact, in order to 
exponentiate H[N] and to expand the exponential in powers, we will have to deal with products of H[Nys. In order 
for these products to be well defined, the domain of the operator must include its range, which is (essentially) Haux- 
Therefore we need to extend the action of H[N] from Hdiff to Haux- The price for this extension is, of course, that 
the operator becomes dependent on the regulator, namely on the precise position in which Cadded is added. However, 
we can do so /lere, because such dependence will disappear in the integration over N ! 

We define H[N] on Haux by simply picking a particular position for Cadded in the definition of Sf 

H[N]\S)^N{x,)A'^{S)\S:). (40) 

where, clearly, A"{S) = Af{s{S)). (When a = (/3, — 1), no modification is necessary. That is, a link is removed 
irrespectively from its precise location.) In a quadratic expression this yields 

{S\H[N]H[N]\s) ^ A^St:) AnS) iV(xs^i=.,J N{xs'^^.,,) (41) 

•1 »2 ■ n 

Here ii labels the nodes of s. After the action of the first operator, and thus the addition (or subtraction) of one link, 
we obtain s"^ . The index 12 runs over the nodes of s"^^ , and therefore its range is larger (or smaller) than the index 
ii, because of the two new nodes (or the two nodes removed) . 



Notice that in each of the terms of the sum in the r.h.s. of (41) (that is, for each fixed value of the indices ii, ai, 12, q;2) 
we have a product 

N{y) N{x) (42) 

where x and y are the positions of the two nodes acted upon by the two operators. In particular, the second H[N] 
operator may act on one of the nodes created by the first operator H[N]. For instance, x in ( ^2[ ) may be the coordinates 
of the point i in (the l.h.s. of) Figure (1), and y may be the coordinates of the point i' (in the r.h.s. of the Figure). 
Now, later on, expressions such as ( ^ ) will appear within functional integrals over N . Inside these integrals, the only 



feature of these two positions that matter is whether x — y or not (see section III). Therefore, the only feature of the 
position of Cadded that matters is whether its end points, namely i' and i" in Figure 1, are on top of i or not. This 
is the only dependence on the regulator (the position of Cadded) that survives in the integral ! More precisely, in the 
integration over N, the arbitrariness in the regularization reduced to the arbitrariness of the decision of whether or 
not we should think at i, i' and i" in the r.h.s. of Figure 1, as on top of the point i on the l.h.s. or not. 

We can view this choice in the following terms. The operator C[N] creates new nodes at positions which are 
displaced from the original node by a distance e, where e is to be later taken to zero (taking this limit is in fact 
necessary in order to identify the quantum operator with the desired classical quantity). The choice is whether to 
take e to zero before or after the integration over N. 

Let us denote the position of the node i of S* (the node acted upon) by x. Denote the two new nodes created by 
the action of the operator as y' and y" . And denote the position of the node i after the action of the operator as y 
(nothing forces x — y a. priori). The natural choices are 

'^^ y = x,y' ^ X, y" X, 

2- 2/ ^ 2:, ?/V X, y" X, 
3. y = X, y' = X, y" = x. 

The choice is exquisitely quantum field theoretical: we are defining here the product of operator valued distributions, 
and we encounter an ambiguity in the renormalization of the regularized product. We thus have three options for the 
regularization of the operator products C[N] . . . C[N], corresponding to the three choices above. 

Choice 3 is not (easily) compatible with the symmetrization of the operator, and choice 2 yields a nonsensical 
vanishing of all the matrix elements of the projector. Thus, we adopt, at least provisionally, choice 1 (which, after 
all, is probably the most natural). That is, we assume that i itself is not displaced by the hamiltonian constraint 
operator, while i' and i" are created in positions which are distinct from the position of i (see Figure 1). 

For a product of n operators (with the same smearing function), we have 

{s\iH[N]r\s)^N{x^)...N{x^) (5) . . . (C::C:r ) (43) 

where ij runs over the nodes of s'"^^ and we have denoted simply the positions of the sequence of 

nodes acted upon in a given term. According to the regularization chosen, this sequence contains points which are 
distinct except when a node is acted upon repeatedly. 
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C. Expansion 



Our task is now to define the space Hphys, using the various tools develo ped above. We aim at defining Hphys 
following the lines Hdiff is defined by the operator 11 in equations (19) and (|l^). That is, we want to construct the 
operator 

P = jDNe~'"^^\ (44) 

whose matrix elements 



{s\P\s') = DN (s|e-'-f^[^l|s'). (45) 



define the quadratic form 



{sW)phys = {S\PW). (46) 



Hphys is th en th e Hilbert space defined over the pre-Hilbert space Hdiff by the quadratic form ( | )phys- As for 11 
(see section p B| ), we will call P a "projector", slightly forcing the usual mathematical meaning of this term. 

Notice that the Hamiltonian constraint we use is a density of weight one (instead than two, as in the original 
Ashtekar formalism); therefore the integration variable iV is a scalar field. This fact will allow us to interpret the 
integration in N in terms of the integral defined in section III J The importance of having a weight-one hamiltonian 
constraint in the quantum theory, was realized by Thiemann H . 

We begin by regularizing the integral (^|) by restricting the integration domain of the functional integral in [DN] 
to the subdomain formed by all the functions N that satisfy 

\N{x)\<T (47) 

where T is a regularization parameter with the dimensions of a time. The physical limit is recovered for T ^ oo. We 
write 

{s\Pt\s')= [ DN {s\e~'"^^^\s'). (48) 

J\N{x)\<T 

Notice that the regularization ( ^7| ) is diffeomorphism invariant. 

By taking advantage from the diff invariance of the expression (p8[), we can insert an integration over the diffeo- 
morphisms and rewrite (^8|) using (^^ as 

{s\Pt\s')^ f D^ f DN {Ui(f>)S\e-'"^^^\s') (49) 

Joiff J\N(x)\<T 



where S is any spin network such that 



Next, we expand the exponent in powers 



s{S) = s. (50) 



{s\PtW) ^ [ Deb [ DN m)S\ ( £ L^{H[N]r] \s'). (51) 

Jozff J\N{x)\<T \„^o n\ J 

Using the explicit form (Esl) of the hamiltonian constraint operator and acting with W((/>) explicitly we obtain 



^One might be puzzled by the fact that the measure defined in section III is normalized, while the measure in equation (pi[), 
which is the formal analog of the expression (ji^), must not be normalized, nor can be seen as the limit of normalized measures. 
The problem, however, is that the choice of the measure in ( ^^ must incorporate the renormalization of the divergence coming 
(at least) from the volume of the gauge orbit. The normalization of the measure is needed to make our expressions converge, 
and should be viewed, we think, as a quantum field theoretical subtraction. 
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n=0 •'D^ff J\N{x)\<T 

xA^(C:-::r)---<^(^)(^lC::r)> (52) 

where ij runs over the nodes of s'"_^"'"^_'^^ . We have also used (^0|) and 

= (53) 

which follows from it. Notice that, as promised, the only remaining diff-dependent quantities are the arguments of the 
functions N{ ). But since the DN integral is diff invariant (see Eq. (|30|)), the integration over Diff can be trivially 
performed using (fsl). Also, notice that the iV(x)'s appear only in the polynomials. Thus we have 



{s\pt\s') = y. A^(C::C:r)•■•<^^) (54) 



71! 

n=0 

where 

4.„....,JT) = / DN N{x,J...N{x,,). (55) 

J\N{x)\<T 

Now, the last integral is precisely the integral of a polynomial Gel'fand function discussed in the previous section. 
The only difference here is that the domain of the dN integral is between —T and T instead than between and T. 
The effect of this is just to put all the odd terms to zero and to double the even terms. Let Uk be the number of 
points that appear k times in the list Xi^ . . . Xi-^ , so that 



kuk = n. (56) 



We obtain 



(57) 



where dni...nt, is defined in Eq. (E9), and e(fc) is defined for any integer n by 



e(2n) = 2, 

e(2n+l) = 0. (58) 



From (p9|), and (|5^), we have 



k ^ k 

Inserting (p9|), in (^), we conclude 

oo 

(s\Pt\s) = Y^ (s|P^"^|s') (60) 

ii . . .ijiOti . . .ctn k 

We recall that the technique for the explicit computation of the coefficients Af{s) is given in |3^]. The last equation 
is an explicit and computable expression, term by term finite, for the regularized matrix elements of the projector on 
the physical state space of the solutions of the hamiltonian constraint. 
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D. Interpretation: spin foam 



The terms of the sum (|6^) are naturally labeled by branched colored surfaces |^,|T^,^, or "spin foams". Each 
surface represents a history of the s-knot state. More precisely, consider a finite sequence ct„ of n + 1 spin networks 

So,Si...,Sn (62) 

In particular, let the sequence ( |6^ ) be generated by a sequence of n actions of single terms of the Hamiltonian 
constraint acting on sq 

'^n^is, <s (63) 

We call such a sequence a "spin foam", and we represent it as a branched colored 2d surface. A branched colored 
surface is a collection of elementary surfaces (faces) carrying a color. The faces join in edges carrying an intertwiner. 
The edges, in turn, join in vertices. A branched colored surface with n vertices can be identified with a sequence ( |6^ ) 
if it can be sliced (in "constant time" slices) such that any slice that does not cut a vertex is one of the spin networks 
in (|63|). In other words, the branched colored surface can be seen as the spacetime world-sheet, or world- history of 
the spin network that evolves under n actions of the hamiltonian constraint. 

Each action of the hamiltonian constraint splits a node of the spin network into three nodes (or combine three nodes 
into one), and thus generates a vertex of the branched surface. Thus, as in the usual Feynman diagrams, the vertices 
describe the elementary interactions of the theory. In particular, here one sees that the complicated action of the 
hamiltonian displayed in Figure 1, which makes a node split into three nodes, corresponds to the simplest geometric 
vertex. Figure 2 is a picture of the elementary vertex. Notice that it represents nothing but the spacetime evolution 
of the elementary action of the hamiltonian constraint, given in Figure 1. 




FIG. 2. The elementary vertex. 



An example of a surface in the sum is given in Figure 3. 
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FIG. 3. A term of second order. 



We write dcjn = s U s' to indicate that the spin foam (t„ is bounded by the initial and final spin networks s and s' . 
We associate to each ct„ the amplitude 

v=l k ^ ' 

where v run over the vertices of (t„ and the amplitude of a single vertex is 

A{v)^j^:{^^:::^). (65) 

The amplitude of a vertex depends only on the coloring of the faces and edges adjacent to the vertex. 
Using this, we can rewrite Eqs. ( |60| , |6l| ) as 

oo 

{s\Pt\s') = Y^ T" ^K)- (66) 

The key novelty with respect to Q is the factor 

k 

The integers Uk are determined by the number of multiple actions of H[N\ on the same vertex. 

The last expression leads immediately to the form of the (regularized) "vacuum to vacuum" transition amplitude, 
or the partition function of the theory 

oo 

n— (7n 

for ct's with no boundaries. In words, the theory is defined as a sum over spin foams sigma„^ where the amplitude 
A((T„) of a spin foam is determined, via (pi), by the product of the amplitudes A{v) of its vertices. Thus, the theory 
is determined by giving the amplitude A{v) of the vertex, as a function of adjacent colors. 
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E. Physical observables 



The expressions we have defined depend on the regulator T . A naive hmit T — > oo yields to meaningless divergences. 
On the other hand, it is natural to expect to be able to remove the regulator only within expressions for physical 
expectation values. Loosely speaking, the integral (^4|) defines a delta-like distribution, and does not converge to any 
function itself; however, its contraction with a smooth function should converge. In particular, for finite T the integral 
contracted with a function gives the integral of the Fourier transform of the function over the interval [— T, T]. If the 
function has a Fourier transform that decays reasonably fast, then the integral should converges nicely. Thus, we may 
expect the expansion in T to be meaningful for suitable observables (see next section). 

The difficulty of constructing interesting physical observables invariant under four-dimensional diffeomorphisms 
in general relativity in well known p3|] and we do not discuss this problem here. Instead, we notice that given an 
operator A on Hdifj, invariant under three-dimensional diffeomorphisms, one can immediately construct a fully gauge 
invariant operator O simply by 

= P A P. (69) 
For instance, A may be the volume of S operator or the projector on a given eigenspace of V 

A = 5{V, v) (70) 
where v is one of the eigenvalues of V. Consider the expectation value of O in a physical state 

, _ {S\0\s)phys _ {S\PAP\S) 

While we expect this quantity to be finite (for an appropriate A) , the numerator and the denominator are presumably 
independently divergent, as one may expect in a field theory. Our strategy to compute (O), therefore, must be to take 
the T — > cxD limit of the ratio, and not of the numerator and of the denominator independently. We thus properly 
define 



and 



(O) - hm {Ot). (73) 

1 — *oo 



Both the numerator and the numerator in ( [7^ ) can be written as power series in T . Therefore we have 

(O) = lini (74) 

where 

an= E {AV\s") (s"|P^"-")|s') (75) 

m— 0,n 

and 

&«-(s|Pj,"^|s). (76) 

The matrix elements {s\P^\s) are explicitly given in (|6l|). Notice that they are finite and explicitly computable. 
Equation ( [72| ) defines a function of T analytic in the origin. We leave the problem of determining the conditions 
under which the higher order terms are small, and of finding techniques for analytically continuing it to infinity on 
the Riemann sphere, for future investigations. 
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F. Quantum ADM surfaces 



An important lesson is obtained by writing the expression for the expectation values in the spin foam version. 
Consider, for simplicity, the case in which A is diagonal in the loop basis. (This is true for the volume, which is the 
reason for the choice of the basis in the intertwiners space, in section |l|). In this case, ( |75| ) becomes 

a„= 5] (sl^'")!^'} A(s')(.s'|P^"-™'|s). (77) 

m— 0,n 

Recalling (^) and this can be rewritten as 

«"=E(E^(^')j^('^n) (78) 

where the s' are all possible spin networks that cut the surface <t„ into two (past and future) parts. 

Equation (^8|) shows that the expectation value of O is given by the average of A on all the discrete ADM-like 
spatial slices s' that cut the quantum spin foam. Summarizing, 

^ ,„ E„r-E.JE.,A(.-))^K) 

The sum in (t„ is over spin foams (with n vertices). For every spin foam, the sum in s' is over all its "spacelike" slices. 

This is a nice geometric result. It clarifies the physical interpretation of the four-dimensional space generated 
by the expansion in T: it is the quantum version of the four-dimensional spacetime of the classical theory. To see 
this, consider for instance the observable defined in (|7^), namely the projector on a given eigenspace of the volume. 
Classically, the volume is defined if a gauge-fixing that identifies a spacelike "ADM" surface S is given. In (^9|) , we see 
that in the quantum theory the role of this spacelike surface is taken by the "ADM-like" spatial slices of the quantum 
spin foam s' . Thus, we must identify the surfaces s' on the spin foam with the classical ADM-surfaces (both are gauge 
constructs!), and therefore we must identify the spin foam itself as the (quantum version of the) four-dimensional 
spacetime of the classical theory.^ 

If the spin foam represents spacetime, the expansion parameter T -introduced above simply as a mathematical trick 
for representing the delta function as the integral of an exponential- can be identified as a genuine time variable (it 
has the right dimensions). This fact provides us with an intuitive grasping on the regime of validity of the expansion 
itself. It is natural to expect that (|7^) might converge for observables that are sufficiently "localized in time" . These 
are precisely the relevant observables in the classical theory as well. We illustrate them in the following section. 



G. 4d diff-invariant observables can be localized in time 



Claims that diffeomorphism invariant observables cannot be localized in time can be found in the literature, and 
have generated much confusion. These claims are mistaken. Let us illustrate why a physical general relativistic 
measurement localized in time is nevertheless represented by a diffeomorphism invariant quantity. 

Consider a state of the solar system. The state can be given by giving positions and velocities of the planets and the 
value of the gravitational field on a certain initial ADM-surface - or, equivalently, at a certain coordinate time t. We 
can ask the following question: "How high will Venus be on the horizon, seen from Alexandria, Egipt, on sunrise of 
Ptolemy's 40th birthday?" In principle, this quantity can be computed as follows. First, solve the Einstein equations 
by evolving the initial data in the coordinate time t. This can be done using an arbitrary time-coordinate choice, and 
provides the Venus horizon height h{t). Then, search on the solution for the coordinate time tpt corresponding to the 
physical event used to specify the time (sunrise time of Ptolemy's birthday, in the example). The desired number is 
finally H = h{tpt), which is a genuine diff-invariant observable, independent from the coordinate time t used. The 
quantity H is coordinate-time independent, but it is also well localized in time. 

In practice, there is no need for computing the solution of the equations of motion for all times t. It is sufficient to 
evolve just from t to tpt, and if t and tpt are sufficiently close, an expansion in (tpt — t) can be effective. 



I thank Mike Reisenberger for this observation. 
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The same should happen in the quantum theory. We characterize the state of the solar system by means of the 
(non-gauge- invariant) state |s), describing the system at a coordinate time t. If we are interested in the expectation 
value of an observable H at a time tpt and if tpt is sufficiently close to t, we may then expect, on physical grounds, 
that the expansion (|7^) be well behaved. In other words, we do not need to evolve the spin network state \s) forever, 
if what we want to know is something that happens shortly after the moment in which the quantum state is |s). 
Nevertheless, what we are computing is 4-d diff invariant quantity. 

As an example that is more likely to be treatable in the quantum theory, consider an observable such as the volume 
Vc of the constant-extrinsic-curvature ADM slice with given extrinsic curvature K(x) = k. Assume we have an 
operator K corresponding to the local extrinsic curvature. Then, 

Vfc - P 5{K, k) V S{K, k) P. (80) 

Once more, Vc is a 4d-diffeomorphism invariant observable localized in time. Given an extrinsic curvature operator 
K, the methods developed here should provide an expansion for the expectation value of Vk- Inserting (|8^) in ([7S|), 
with O = Vk, the delta function selects the ADM slices with the correct extrinsic curvature from the sum in s', and 
the mean value of Vk is given by the average over such slices appearing in the time development of |s) generated by 
the Hamiltonian constraint. If \s) is sufficiently close (in time) to a K = k surface, then, on physical grounds, we 
have some reasons to hope that the expansion T to be well behaved. Thus, the diffeomorphism invariant quantum 
computation of Vk reproduces the structure of the classical computation. 



V. CONCLUSIONS 



We have studied the dynamics of nonperturbative quantum gravity. Because of the diffeomorphism invariance 
of the theory, this dynamics is captured by the "projector" P on the physical states that solve the hamiltonian 
constraint. We have constructed an expansion for the (regularized) projector P, and for the expectation values of 
physical observables. The expansion is constructed using some formal manipulations and by using a diffeomorphism 
invariant functional integration on a space a scalar functions. This construction may represent a tool for exploring 
the physics defined by various hamiltonian constraints. 

Our main result is summarized in equations ( ^0| - |6l| ), which give the regularized projector and equations (|7^[76|), 
which gives the expectation value of a physical observable, both in terms of finite and explicitly computable quantities. 
Equivalently, the theory is defined in the spin-foam formalism by the partition function Z given in equations ( [64[ - |65| - 
|68| ) . The expectation values are then given in equation ( |79| ) as averages over the spin foam. The spin- foam formalism 
is particularly interesting, because it provides a spacetime covariant formulation of a diffeomorphism-invariant theory. 
The partition function Z is expressed "a la Feynman" as a sum over paths, but these paths arc topologically distinct, 
and discrete (so that we have a sum rather than an integral). 

Several aspects of our construction are incomplete and deserve more detailed investigations, (i) The physical 
discussion on the range of validity of the expansions considered should ce rtain ly be made more mathematically 



precise, (ii) The choice of the position of the nodes considered in Section IV B is somewhat arbitrary and other 
possibilities might be explored, (iii) A possible modification of the formalism that could be explored is to restrict the 
range of the integration in N to positive definite A'^'s, in analogy with the Feynman propagator, (iv) The spacetime 
geometry of the (individual) spin foams has not yet been fully understood and deserves extensive investigations (on 
this, see ||l7| , [l8| , ^ , p2| , ^ ) . (v) We have completely disregarded the Lorentzian aspects of the theory. These can be 
taken into account by using the Barbero-Thiemann Lorentzian hamiltonian constraint Q, or, alternatively, along the 
lines explored by Smolin and Markopoulou 0. (vi) The limit T — > oo should be better understood: how can we find 
it from the knowledge of a finite number of the a„ and 6„ coefficients in (^)? (vii) Finally, the formalism developed 
here makes contact with the spin foam models pO|-p^ ) We believe that the relation between these approaches deserves 
to be studied in detail. 

The key issue is whether the measure that we have emploicd is the "correct" one. Intuitively, whether this measure 
has the property that the integral of the exponential gives the delta function, or whether P is in fact, in the appropriate 
sense, a projector. We will discuss this point elsewhere. 

The nonperturbative dynamics of a diffeomorphism invariant quantum field theory is still a very little explored 
territory; the scheme proposed here might provide a path into this unfamiliar terrain. 

I thank Don Marolf, Mike Reisenberger, Roberto DePietri, Thomas Thiemann and Andrea Barbieri for important 
exchanges and for numerous essential clarifications. This work was supported by NSF Grant PHY-95-15506, and by 
the Unite Propre de Recherche du CNRS 7061. 



15 



[1] 

[2] 
[3] 

[4] 

[5] 
[6] 
[7] 
[8] 



[9] 
[10] 

[11] 



[12] 
[13] 
[14] 

[15] 

[16] 
[17] 
[18] 
[19] 
[20] 
[21] 

[22] 

[23] 



[24] 
[25] 

[26] 

[27] 
[28] 
[29] 

[30] 
[31] 
[32] 
[33] 
[34] 

[35] 

[36] 



C Rovelli, L Smolin Phys. Rev. Lett. 61, 1155, (1988); Nucl. Phys. B331 (1), 80-152, (1990). 

A Ashtekar, Phys. Rev. Lett. 57 (18), 2244-2247, (1986); Phys. Rev. D36 (6), 1587-1602, (1987). 

For an introdu ction and complete references, see C Rovelli "Loop Quantum Gravity " , Living Revie ws in Relativity, elec- 
tronic journal, http://www.livingreviews.org /Articles / Volumel / 1998- Irovellj ; also in ^r-qc/97090"08l 

C Rovelli, L Smolin, Nucl. Phys. B442, 593-622, (1995). Erratum: Nucl. Phys. B456, 734, (1995). A Ashtekar, J 
Lewandowski, Class, and Quantum Grav. 14, A55-A81, (1997). 

T Thiemann, Phys. Lett. B380, 257-264, (1996); Class, and Quantum Grav. 15, 839 (1998) . 
M Reisenberger, C RoveUi, Phys. Rev. D56, 3490-3508 (1997). 

F Markopoulou, L Smolin, Nucl. Phys. B508, 409-430 (1997); Phys Rev D58, 84033 (1998). 

J Lewandowski, D Marolf, Int. J. Mod. Phys. D7, 299-330 (1998). R Gambini, J Lewandowski, D Marolf, J PuUin, 

Int. J. Mod. Phys. D7, 97-109, ( 1998). L Smoli n, "The classical limit and the form of the hamiltonian constraint in 

nonperturbative quantum gravity", j;r-qc/9609034. 

C RoveUi, J. Math. Phys. 36, 5629 (1995). 

J Klauder, Annals Phys. 254, 419-453 (1997). 

J Govaerts, "Projection Operator Approach to Constrained Systems" J. Phys. A 30, 603-617 (1997). J Govaerts, J 
Klauder, "Solvin g Gauge Invariant Systems without Gauge Fixing: the Physical Projector in 0-1-1 Dimensional Theories", 
hep-th/9809lT9| . 

LV Prokhorov, SV Shabanov, Phys Lett B216, 341 (1989); Sov. Phys. Uspekhi 34, 108 (1991). 

M Henneaux, C Teitelboim, "Quantization of constrained systems" (Princeton University Press, 1992). 

SW Hawking, in Relativity, Groups and Topology, Les Houches Session XL, B DeWitt, R Stora eds, (North Holland, 

Amsterdam, 1984). 

A Ashtekar, J Lewandowski, J. Math. Phys. 36, 2170 (1995). JC Baez, Lett. Math. Phys. 31, (1994); Lett. Math. Phys. 
31, 213-224, (1994). 

T Thiemann, Class. Quantum Grav. 15 (1998) 1487-1512. 
J Baez, Class. Quantum Grav. 15 (1998) 1827-1858. 

M Reisenberger, "Worldsheet formulations of gauge theories and gravity", gr-qc/9412035. 

J Iwasaki, J. Math. Phys. 36, 6288, (1995). 

J Barret, L Crane, J. Math. Phys., 39, 3296-3302 (1998). 

M Reisenberge r, Class. Quantum Grav., 14, 1753-1770; "A lattice worldsheet sum for 4-d Euclidean general relativity", 
gr-qc/971105^ . 



M Reisenberger, "Classical Euclidean general relativity from "left-handed area 
De Pietri, L Freidel, "so(4) Plebanski Action and Relativistic Spin Foam Model", |gr-qc/9804071 



right-handed area", gr-qc/9804061 



R 



For an overview of the present approaches to the problem of the description of quantum spacetime, see: C Rovelli, "Strings, 
loops and the others: a critical survey on the present approaches to quantum gravity", in Gravitation and Relativity: At 
the tu rn of the Millen ium. Proceedings of the GR15, Pune, India 1997, N Dadhich J Narlikar eds (Puna University Press, 
1998), |gr-qc/980302^ 

Loop quantum gravity and quanta of space: a primer", gr-qc/9806079[ 

D49, 6935-6938, (1994); Phys. Rev., D 51, 5507-5510, (1995); Phys. Rev., D 51, 5498-5506, 



Class, and Quantum Grav., 9, 1433-85, (1992). A Ashtekar, J Lewandowski, D Marolf, J Mourao, 
Phys. 36, 6456-6493, (1995). 

Phys. Rev. D54, 2664-2690, (1996). R DePiet ri, Class, and Q uantum Grav 14, 53-69, (1997). 



C Rovelli, P Upadhya, 
F Barbero, Phys. Rev 
(1995). 

A Ashtekar, CJ Isham 
T Thiemann, J. Math. 
R DePietri, C Rovelli, 
JC Baez, S Sawin, "Diffeomorphism-Invariant Spin Network States", :j-alg/9708005 . 

CJ Isham, in Relativity, Groups and Topology II, Les Houches 1983, B DeWitt and R Stora eds (North Holland Amsterdam 
1984). 

NP Landsman, C/ass and Quantum Grav 12, L119 (1995). 

R Borissov, R DePietri, C Rovelli, Class and Quantum Grav 14, 2793, (1997). 
C RoveUi, L Smolin, Phys Rev Lett 72, 1994 (446). 

C Rovelli, Class and Quantum Grav 8, 297 (1991); Class and Quantum Grav 8, 317 (1991). 

R Loll, Phy s. Rev. Lett. 7 5, 3048-3051, (1995). A Ashtekar, J Lewandowski, "Quantum Theory of Geometry II: Volume 



operators", gr-qc/9711031 . 

C Rovelli, L Smolin, Phys. Rev. D52, 5743-5759, (1995). JC Baez, Adv. m Math. 117, 253, (1996); in The Interface of 
Knots and Physics, LH Kauffman ed, (American Mathematical Society, Providence, Rhode Island, 1996), pp. 167-203. 
A Barbieri, Nucl. Phys. B518, 714, (1998). 



16 



